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Abstract 

We discuss regularization by noise of the spectrum of large random non- 
Normal matrices. Under suitable conditions, we show that the regularization 
of a sequence of matrices that converges in *-moments to a regular element 
a, by the addition of a polynomially vanishing Gaussian Ginibre matrix, 
forces the empirical measure of eigenvalues to converge to the Brown mea- 
sure of a. 

1 Introduction 

Consider a sequence A N of N xN matrices, of uniformly bounded operator 
norm, and assume that An converges in *-moments toward an element a 
in a W* probability space (fL,\\ • ||,*,cp), that is, for any non-commutative 
polynomial P, 

^trP{A N ,A* N ) q>(P(a,a*)) . 

We assume throughout that the tracial state cp is faithful; this does not rep- 
resent a loss of generality. If An is a sequence of Hermitian matrices, this 
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is enough in order to conclude that the empirical measure of eigenvalues of 
An, that is the measure 



jA ._ 



1 N 

-Yh. 



where X,(A^),/ = I...N are the eigenvalues of An, converges weakly to 
a limiting measure /u a , the spectral measure of a, supported on a compact 
subset of E. (See 0] Corollary 5.2.16, Lemma 5.2.19] for this standard 
result and further background.) Significantly, in the Hermitian case, this 
convergence is stable under small bounded perturbations: with Bn = An + 
En and ||£V|| < e, any subsequential limit of L N will belong to B^^, 8(e)), 
with 8(e) — ^e^o and Br,{\>a,r) is the ball (in say, the Levy metric) centered 
at v fl and of radius r. 

Both these statements fail when A„ is not self adjoint. For a standard 
example (described in @), consider the nilpotent matrix 



in 
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Obviously, L N = So, while a simple computation reveals that Tn converges 
in *-moments to a Unitary Haar element of A , that is 



tr(T^(T*)^...T^(T*)^) 



N 



0, otherwise. 



(1) 



Further, adding to Tn the matrix whose entries are all except for the bot- 
tom left, which is taken as e, changes the empirical measure of eigenvalues 
drastically - as we will see below, as ,/V increases, the empirical measure 
converges to the uniform measure on the unit circle in the complex plane. 

Our goal in this note is to explore this phenomenun in the context of 
small random perturbations of matrices. We recall some notions. For a G A , 
the Brown measure v a on C is the measure satisfying 

logdet(z-a) = j \og\z- z\dv a (z), z e C, 

where det is the Fuglede-Kadison determinant; we refer to O SJ for defini- 
tions. We have in particular that 



log det (z — a) 



\ogxdv 7 a (x) z€C. 
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where v~ denotes the spectral measure of the operator \z — a\. In the sense 
of distributions, we have 

v a = -!-Alogdet(z-£?). 
271 

That is, for smooth compactly supported function \|/ on C, 

J \\r(z)dv a (z) 



A crucial assumption in our analysis is the following. 
Definition 1 (Regular elements). An element a G A is regular if 

lim / dzA\\f(z) / logxdv^x) = 0, (2) 
e ^°Jc Jo 

for all smooth functions \|/ on C with compact support. 

Note that regularity is a property of a, not merely of its Brown measure 
v a . We next introduce the class of Gaussian perturbations we consider. 

Definition 2 (Polynomially vanishing Gaussian matrices). A sequence of 
N-by-N random Gaussian matrices is called polynomially vanishing if its 
entries (Gn(i, j)) are independent centered complex Gaussian variables, and 
there exist K > 0, k' > 1 + K so that 

AT* <E\Gij\ 2 <N^ K . 

Remark 3. As will be clear below, see the beginning of the proof of Lemma 
[lOl the Gaussian assumption only intervenes in obtaining a uniform lower 
bound on singular values of certain random matrices. As pointed out to us by 
R. Vershynin, this uniform estimate extends to other situations, most notably 
to the polynomial rescale of matrices whose entries are i.i.d. and possess a 
bounded density. We do not discuss such extensions here. 

Our first result is a stability, with respect to polynomially vanishing 
Gaussian perturbations, of the convergence of spectral measures for non- 
normal matrices. Throughout, we denote by [|M [| p the operator norm of a 
matrix M. 



— J dz&y{z) j \og\z-z\dv a {z) 

— j dzAy(z) J logxdv z a (x). 
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Theorem 4. Assume that the uniformly bounded (in the operator norm) se- 
quence ofN-by-N matrices An converges in ^-moments to a regular element 
a. Assume further that L N converges weakly to the Brown measure \ a . Let 
Gm be a sequence of polynomially vanishing Gaussian matrices, and set 
B N = A N + Gn- Then, L N — > V a weakly, in probability. 

Theorem [4] puts rather stringent assumptions on the sequence An- In 
particular, its assumptions are not satisfied by the sequence of nilpotent ma- 
trices Tn in {J}- Our second result corrects this defficiency, by showing that 
small Gaussian perturbations "regularize" matrices that are close to matrices 
satisfying the assumptions of Theorem [4] 

Theorem 5. Let A^, En be a sequence of bounded (for the operator norm) 
N-by-N matrices, so that An converges in ^-moments to a regular element 
a. Assume that \\E^\\ op converges to zero polynomially fast in N, and that 
L N +E —7- V a weakly. Let Gn be a sequence of polynomially vanishing Gaus- 
sian matrices, and set Bn = An + Gn- Then, L N — > \ a weakly, in probability. 

Theorem [5] should be compared to earlier results of Sniady [6], who 
used stochastic calculus to show that a perturbation by an asymptotically 
vanishing Ginibre Gaussian matrix regularizes arbitrary matrices. Compared 
with his results, we allow for more general Gaussian perturbations (both 
structurally and in terms of the variance) and also show that the Gaussian 
regularization can decay as fast as wished in the polynomial scale. On the 
other hand, we do impose a regularity property on the limit a as well as on 
the sequence of matrices for which we assume that adding a polynomially 
small matrix is enough to obtain convergence to the Brown measure. 

A corollary of our general results is the following. 

Corollary 6. Let Gn be a sequence of polynomially vanishing Gaussian 
matrices and let Tn be as in ([TJ). Then L^ +G converges weakly, in probability, 
toward the uniform measure on the unit circle in C. 

In Figured] we give a simulation of the setup in Corollary [6] for various N. 

We will now define class of matrices Tj, jv for which, if b is chosen cor- 
rectly, adding a polynomially vanishing Gaussian matrix Gn is not sufficient 
to regularize 7& ^ + Gn- Let b be a positive integer, and define 7& ^ to be an 
Nby N block diagonal matrix which each b + 1 by b + 1 block on the diag- 
onal equal 7^+1 (as defined in ([TJ). If b + 1 does not divide Af evenly, a block 
of zeros is inserted at bottom of the diagonal. Thus, every entry of Tt, is 
zero except for entries on the superdiagonal (the superdiagonal is the list of 
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(a) Af = 50 




(b) N = 100 




(c) N = 500 




(d) N = 5000 
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Figure 1 : The eigenvalues of Tm + N 3 
with mean 0, variance 1 entries. 


1//2 Gjy, where Gn is iid complex Gaussian 



entries with coordinates 1) for 1 < i < N — 1), and the superdiagonal 
of r^ jv is equal to 

°' 1 > 1 > • • • ' 1 ' °' • • • ' 1 > 1 > • • ■ > 1 > o^^^o) • 

b b b <b 

Recall that the spectral radius of a matrix is the maximum absolute value of 
the eigenvalues. Also, we will use |A|| = tr(A*A) 1 / 2 to denote the Hilbert- 
Schmidt norm. 

Proposition 7. Let b = b(N) be a sequence of positive integers such that 
b(N) > logN for all N, and let TbN be as defined above. Let R^ be an 
N by N matrix satisfying \\Rn\\ < <?(A0, where for all N we assume that 



>(N)< 



3b^N 



. Then 



p(Tb,N + RN)<{Ng(N)) 1/h + o(\) 



where p(M) denotes the spectral radius of a matrix M, and o(l) denotes a 
small quantity tending to zero as N — )■ °°. 

Note that Ty^ converges in *-moments to a Unitary Haar element of A 
(by a computation similar to (Q])) if b(N) /N goes to zero, which is a regular 
element. The Brown measure of the Unitary Haar element is uniform mea- 
sure on the unit circle; thus, in the case where {Ng(N)) l l b < 1, Proposition |7] 
shows that 7&jv +Rn does not converge to the Brown measure for Tj,^. 

Corollary 8. Let Rfj be an iid Gaussian matrix where each entry has mean 
zero and variance one. Set b = b(N) > logAf be a sequence of integers, and 
lety> 5/2 be a constant. Then, with probability tending to 1 as N —> °°, we 
have 

p(T b , N + exp(-yb)R N ) <exp(-Y+— ^— J +o(l), 
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(a) N = 50 




(b) N = 100 




(c) N = 500 




(d) N = 5000 
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Figure 2: The eigenvaules of T[ og N,N+N 


3 ~ l/2 G N , where G N 


is iid complex 


Gaussian with mean 0, variance 1 entries. The spectral radius is roughly 0.07, and 


the bound from Corollary [8] is exp(— 1) 


^0.37. 









where p denotes the spectral radius and where o( 1) denotes a small quantity 
tending to zero as N — >■ °°. Note in particular that the bound on the spectral 
radius is strictly less than exp(— 1/2) < 1 in the limit as N — > °°, due to the 
assumptions on y and b. 

Corollary [8] follows from Proposition |7] by noting that, with probability 
tending to 1, all entries in are at most ClogiV in absolute value for some 
constant C, and then checking that the hypotheses of Proposition [7] are sat- 
isfied for g(./V) = exp(— yZ?)C/V(log./V) 1 / 4 . There are two instances of Corol- 
lary [8] that are particularly interesting: when b = N — 1, we see that a expo- 
nentially decaying Gaussian perturbation does not regularize = 7]v_i ^, 
and when b = log(Af), we see that polynomially decaying Gaussian pertur- 
bation does not regularize T\ og ^^ (see Figure [2]). 

We will prove Proposition [7j in Section [5] The proof of our main results 
(Theorems[4]and[5]) borrows from the methods of [3 ]. We introduce notation. 
For any N-by-N matrix C#, let 



N 



o c N 
Q o 



We denote by Gq the Cauchy-Stieltjes transform of the spectral measure of 
the matrix Cn, that is 



Gc(z) 



1 

2N 



tr(z-Cv) 



-l 



ZGC+. 



The following estimate is immediate from the definition and the resolvent 
identity: 



\G c (z)-G D (z)\< llC ^ op 



(3) 
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2 Proof of Theorem a 



We keep throughout the notation and assumptions of the theorem. The fol- 
lowing is a crucial simple observation. 

Proposition 9. For all complex number % and all z so that 3z > N~ 5 with 
8 < k/4, 

E|3G BwH ( Z )|<E|3G A ^(z)| + l 

Proof. Noting that 

K\\B N - A N \\ k op = K\\G N \\ k op < CkN-^ 2 , (4) 

the conclusion follows from Q and Holder's inequality. □ 
We continue with the proof of Theorem [4] Let denote the empirical 

measure of the eigenvalues of the matrix An — Z- We have that, for smooth 
test functions \|/, 

J dzAy{z) j \og\x\d\\ N {x) = ^J V(z)dl4f(z). 

In particular, the convergence of L N toward v a implies that 

E J dzA\\r(z) J log \x\dv% N (x)-> J y(z)dv a (z) = J dzA\\r(z) J logxdv z a (x) 

On the other hand, since x i-> log x is bounded continuous on compact subsets 
of (0,o°), it also holds that for any continuous bounded function £ : M + \-t R 
compactly supported in (0,°°), 

E J dzA^{z) j t,{x)\ogxdv\ N (x) -> J dzbty(z) j t,{x)\ogxdv z a {x) . 

Together with the fact that a is regular and that A^ is uniformly bounded, 
one concludes therefore that 



lim lim E / / log\x\dv\ (x)dz = 0. 
40^» J Jo An 

Our next goal is to show that the same applies to B^. In the following, we 
let v z Bn denote the empirical measure of the eigenvalues of Bn — z. 

Lemma 10. 

lim lim / E[ / log Ixl^dvi (x)]dz = 
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Because E||fiw — Awllop ~~ ^ for any k > 0, we have for any fixed smooth 
w compactly supported in (0, °°) that 



?N^°° O , 



E| J dz,Ay(z) J w(x) logxdv\ N (x)- J dzA\\r(z) J w(x) logxdv z BN (x) 
Theorem |4] follows at once from Lemma [TOl 

Proof of lemma [TOj Note first that by ||5] Theorem 3.3] (or its generalization 
in (3l Proposition 16] to the complex case), there exists a constant C so that 
for any z, the smallest singular value o N of B^+zJ satisfies 

P(a z N <x) <C(A^ +K ^ 



with P = 1 or 2 according whether we are in the real or the complex case. 
Therefore, for any £ > 0, uniformly in z 

E[ 



:[J log|x|-^v| w (x)] < Ellogio^)- 1 ^^} 

= c(^+ K '-^ P log(^) + ^ V 1 -c(n 



i / i.i 

x dx 



goes to zero as N goes to infinity as soon as t, > 4 + k'. We fix hereafter 
such a £ and we may and shall restrict the integration from to £. To 
compare the integral for the spectral measure of An and B^, observe that for 
all probability measure P, with P y the Cauchy law with parameter y 

P([a,b})<P*Py([a-n,b + n])+P y ([-n^Y)<P*P y ([a-r\,b + n}) + ^ 

(5) 

whereas for b — a > r\ 

P([a,b])>P*Py([a + n,b-r\])-^. (6) 

Recall that ^ 

P*P y ([a,b}) = j \ZG{x + q)\dx. (7) 

J a 

Set y = N~ K ^ 5 , k" = k/2 and r| = N^'^ 5 . We have, whenever b — a> 4r\, 

rb+r\ 

Ev% N ([a,b]) < / E|3G b „ +z (x + /y)|</x + aH k - k ">/ 5 

Jfl-T| 

< (6- a + 2N- K "/ 5 ) +v z An *P N -, /5 ([a - Ar K / 10 ,fc+Ar K / 10 ]) + AT K / 10 

< (b-a + 2N- K/w )+v z A ([(a-2N- K/w ) + ,(b + 2N- K/lQ )}) + 2N- K/W , 
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where the first inequality is due to ([5]) and (0, the second is due to Proposi- 
tion|9j and the last uses © and (O. Therefore, if b — a = C/V~ K / 10 for some 
fixed C larger than 4, we deduce that there exists a finite constant C which 
only depends on C so that 

Ev| N ([a^])<C'(ft-a)+v^([(a-2^- K / 10 ) + ,(^ + 2^/ 10 )]). 

As a consequence, as we may assume without loss of generality that k' > 
k/10, 



E[f \o % \x\~ x dN\ N {x)\ 



[jV*/ 10 e] 

< + 2CkN- K / l0 )- l E[v z BN ]([N^ + 2CkN- K / l0 ,N^ + 2C(k + 1)AT K / 10 ]) . 

k=0 

We need to pay special attention to the first term that we bound by noticing 
that 

log(A^)-i E [v|J[AK,AK + 2CAr K / 10 ])] 

< l^log(A^- K / 10 )- 1 E[v| w ([0,2(C+ l)Ar K / 10 ])] 

< -^log(^- K / 10 )- 1 (2C , Af- K / 10 +v^([0,(C + 2)A^- K / 10 ])) 

1C 



20C% 



f 2(C+2)N- K /'° 

log{N- K ^°)- l N- K ^ + C" / tog|*|- 1 dv^(x) 

J o 



For the other terms, we have 

[ATK/10 E ] 

log(^-^ + 2C^ K/10 )- 1 E[v^]([iV-^ + 2C^ K/10 ,^ ? + 2C(^+l)^ K/10 ]) 

k=i 

[yVK/io e ] 

< 2C' iogiCkN-^ 10 )-^^^ 10 
k=\ 

+ Y log(CkN- K/w y i v z AN ([2C{k-l)N- K/m ,2C(k + 2)N- K/l0 }). 
k=\ 

Finally, we can sum up all these inequalities to find that there exists a finite 
constant C'" so that 

f'E f'E f£ 

E [ / log | x\ ~ 1 dv z B (x)] < C" / log \x\ - 1 dv\ (x) + C" / log \x\- 1 ^ 



and therefore goes to zero when n and then £ goes to zero. This proves the 
claim. □ 



3 Proof of Theorem S 

From the assumptions, it is clear that (An + En) converges in *-moments 
to the regular element a. By Theorem 01 it follows that L N +E+G converges 
(weakly, in probability) towards v fl . We can now remove En- Indeed, by © 
and ©, we have for any % < k' /2 and all £, G C 

\ G A+G+^ Z ) ~ G A+G+E+t,(z)\ ^ 

and therefore for 3z > N^ 2 , 

\ZG N A+GH (z)\<\ZG N A+G+E ^(z)\ + l. 

Again by [5, Theorem 3.3] (or its generalization in [3, Proposition 16]) to 
the complex case), for any z, the smallest singular value o z N of An + Gn + z 
satisfies 



p{o N <x) <c[m +K: x 



with P = 1 or 2 according whether we are in the real or the complex case. We 
can now rerun the proof of Theorem |U replacing An by A' N = An +En + Gn 
and Bn by A' N — En- □ 



4 Proof of Corollary [6] 

We apply Theorem [5] with An = Tn, En the N-by-N matrix with 

F (i , w 8 n = n-^\ i = ij=N 

0, otherwise, 

where k' > K. We check the assumptions of Theorem [5] We take a to be 
a Unitary Haar element in A , and recall that its Brown measure v a is the 
uniform measure on {z 6 C : \z\ = 1}- We now check that a is regular. 
Indeed, f x k dv z a (x) = if k is odd by symmetry while for k even, 



/ 



x^ ( (x)=cp([(z-a)(z-a)f/ 2 ) = ^(|z| 2 + l)^( \. ) ( J 
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and one therefore verifies that for k even, 



J x k dv z a (x) = ^ j (|z| 2 + l+2|z|cos6)^e. 



It follows that 

\ogxdv z a (x) = — log(|z| 2 + l+2|z|cose)l { | z |2 +1+2 | z | cos6<e} t/e^ E ^ 0, 

proving the required regularity. 

Further, we claim that L^ +E converges to \ u . Indeed the eigenvalues X 
of An + En are such that there exists a non- vanishing vector u so that 

that is 

In particular, all the Af -roots of 8jv are (distinct) eigenvalues, that is the eigen- 
values M of An are 

X N j =\h N \ l l N e 2i ^l N , l<j<N. 
Therefore, for any bounded continuous g function on C, 



1=1 

as claimed. □ 



5 Proof of Proposition 13 

In this section we will prove the following proposition: 

Proposition 11. Letb = b(N) be a sequence of positive integers, and let Tb^ 
be as in Proposition^ Let Rfj be anNbyN matrix satisfying \\Rn\\ < g(N), 
where for all N we assume that g(N) < 3fc ^y - Then 

P (r^+^)<(o(Vw(2^ 
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Proposition [7]follows from Proposition [TT]by adding the assumption that 
b(N) > log(N) and then simplifying the upper bound on the spectral radius. 

Proof of Proposition [TTJ To bound the spectral radius, we will use the fact 

that p(Tb,N + Rn) < ||(7/?,A' + ^A') <: || 1 ^ f° r all integers k> 1. Our general 
plan will be to bound + -Rat) an d then take a k-th root of the bound. 

We will take k = b + 1, which allows us to take advantage of the fact that 
Tbjf is (b + l)-step nilpotent. In particular, we make use of the fact that for 
any positive integer a, 



I T a I 
\ 1 b.N\ 



N I 1/2 



if 1 < a < b 
if b + l <a. 



(8) 



We may write 

\\(Tb.N' 



R? 



b+l 



i=l 



s E 

%€{0,l} b + l 
b+l 

=E E 

(=o \e{0,i} b + l 

X has I ones 



b+l 



2V 



!=1 



When ^ is large, we will make use of the following lemma. 
Lemma 12. //X e {0,1}^ /ias £ ones and £ > (jfc + l)/2, f/zera 



II ^ 



N 



i=\ 



< 



-!_*-<+ 1 J 



l b,N 



k-l+1 



\Ra 



ik-l 



We will prove Lemma [T2l in Section 157X1 

Using LemmafT2lwith k = b + 1 along with the fact that ||AB|| < ||A|| \\B\\, 
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we have 



(7&,Af+^w) fo+1 || < 



f=0 v 7 



\b-i+\ 



b+\ 



+ E 



&+i\ 



7^ L 5=1+3. 
i 6JV 



li>-M-l 



< E (tMiw' 

£=0 



(9) 



+ 



e CD 



' U-*+2_l 



L b,N 



b-l+2 



8 b - l+ \ (10) 



where the second inequality comes from the assumption \\Rn\\ < g = g(-^0- 
We will bound (O and (fTOl) separately. To bound © note that 



•+1 



\\n, 



l fe+2 

1=0 

b+Af b+l 



+ n ((»+!) 



_/V_ 
6 + 1 



f/2 



J>-^+l 



< 



2 U(* + 1)/2J 
(VNb(2N 1 / 4 g 1 / 2 ) b 



N (b+2)/A b/2 



(11) 



Next, we turn to bounding (ITQb . We will use the following lemma to 
show that the largest term in the sum (fTOl) comes from the l = b term. Note 
that when £ = b+l, the summand in (fTOl) is equal to zero by ([8]). 



Lemma 13. ; If 



rL*-«+l. 



l b,N 



> a«<i t < b — 1 a?i<i 



e V2 N l/2 b > 



then 



,[5=1+2. 



L b,N 



b-M 



< 



V+lJ 

We will prove Lemma [T3l in Section [5+1 



• I. 5=1+1. 



l b,N 



b-£+l 



J>-1 
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Using Lemma [T3l we have 



b+\ 

E 



'+1 



rr< L b-I+2_ 

£ b,N 



b-l+2 



11 



<-( b+ i)(b- 

< b 2 Ng. 



N 



(12) 



Combining CO} and C[2]> with © and (fTOl) . we may use the fact that 
(x + yyK b+ ^ < jcV(^I) +yi/(*+l) for positive x,j to complete the proof of 
Proposition [TT1 It remains to prove Lemma [121 and Lemma [T3l which we do 
in Section [5TT1 below. □ 



5. 1 Proofs of Lemma [12] and Lemma [13] 



Proof of Lemma [I2j Using dHJ, it is easy to show that 



\T a \\\\T C \\ <r 
\ 1 b.N\\ \\ 1 b,N\\ < - 



pa— 1 

b,N 



T c+1 
1 b.N 



for integers 3 < c + 2 < a < b. (13) 



It is also clear from (HJ that 



It" II < 

\ 1 b.N\\ - 



pa— 1 
l b.N 



for all positive integers a. 



(14) 



Let X G {0, 1} have ^ ones. Then, using the assumption that £ > k — t + 
1, we may write 



X, pl-X/ -p«i T>b\rpa2 nb2 rpak-LTfbk-Lrpak-^x 

N b.N N b,N N ' ' ' b,N N 2 b,N ' 



i=l 



where a, > 1 for all i and bj > for all i. Thus 



x, 



/'=1 



n to 



(=i 



Applying (TT31 repeatedly, we may assume that two of the a, differ by more 
than 1, all without changing the fact that Xw=f +1 fl i = Thus, some of the 
at are equal to J and some are equal to Finally, applying 

(fT4l . we have that 



lc-£+\ 

n Ita 



< 



-p Lw+i. 
*b,N 



k-e+i 



14 



□ 



Proof of Lemma [I3j Using © and rearranging, it is sufficient to show that 



£+\ 



b 



b-£ + 2 



+ 1 



1/2 



N 
b+l 



1/2 



< 



Ms* 



e+i. 



+ 1 



+ i 



Using a variety of manipulations, it is possible to show that 

b-e+i 

+ M 2 / (b-£ + 2){b-£+\) 



.6-^+2. 



+ 1 



> 



exp 



& + 2 



(6 + 2)(6-^ + 2)-^ (6 + 2)(6-£ + 2) 
> exp(-3/2). 



Thus, it is sufficient to have 



^ 1 /2 g < eX p(-3/2), 



which is true by assumption. 



□ 
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